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Abstract 

If Ci is the convex hull of the curve of the standard Brownian motion in the 
complex plane watched from to 1, we consider the convex hulls of Ci and sev- 
eral rotations of it and we compute the mean of the length of their perimeter by 
elementary calculations. 

I Introduction 

If C is a compact plane convex set, denote by L{C) the length of its perimeter. Let 
B(t) = {X(t), Y(t)) be a standard Brov^^nian motion in the plane and denote by Ci(t) the 
convex hull of the curve 

B = {B{s) ■,0<s<t}. 

This object has been considered by Paul Levy (1948), pages 239-240. The first author 
a long time ago has called attention on the fact that the mean of L(Ci{t)) is y/Snt in 
Letac (1978). Here we consider the symmetric convex hull 6*2(^)5 namely the convex hull 
of Ci{t) U — Ci(t). The convex set C2(t) should not be confused with the Minkowski sum 
Ci(t) + (— Ci(t)), which is bigger. By scaling we can guess the existence of a constant 
£2 such that E[L{C2{t))] = £2^87?^. Similarly we consider random convex sets obtained 
in the following way. We let fl C [0,27r[ be a set of angles and the convex hull of 
Ucjen-R^B, where i?^ is the rotation of angle u. Again by scaling there exists a constant 
ifi such that E[L{CQ{t))] = ifiy/Sirt. We shall compute these constants when f2 is equal 
to one of the following sets 

{0}, {0,7r}, {0,7r/2}, {0,27r/3}, {0,7r/2,7r} 

or 

{0,27r/3,47r/3} {0, 7r/2, vr, 37r/2}, [0,2tt[ 

giving the result either in closed form or as an expression which can be evaluated numer- 
ically. 
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II Computations 



II. 1 Preliminary remarks 

Given a compact convex set C in the Euclidean plane consider its support function 
he as the function on R of period 2ti defined by 

^c(^) = niax (xcos6' + 'usin6'). 

It is well known (see Letac (1983), formula (7.2) for a proof without regularity hypothesis) 
that the length of the perimeter of C is hc{0)d9. This formula is generally attributed 
to Cauchy. Let Q C [0, 2n[ and define h-B{0) = maxo<s<t(X(s) cos6' + Y{s) sin 6*) and 

hniO) = suphsiO + uj), 

then rotation invariance of the law of Brownian motion implies 

E[L{Cnm = 27rE[hnm- 

Let Pn be the convex set defined by the inequalities 

xcosu + y sin uj < 1; u & Q. 

In other terms, Pq is the intersection of half planes whose boundary is tangent to the 
unit circle. Denote by the first time where Brownian motion exits this set. Then a 
scaling argument shows that /^^(O) is distributed as . It turns out that in several 
cases the Laplace transform i?[exp(— ^T^)] is given by a simple formula. Under such a 
circumstance _ 

- E[exp(--T^)]dA 

E[eM-^Ta)]d\. (2.1) 



and 



II. 2 Some easy cases 

II.2.1 fi = {0} 

The time is the first hitting time of 1 for a linear Brownian motion, and 

£;[exp(-yTo)] = e-\ (2.2) 
The Laplace tranform can be inverted to give 

P{Tn e dt) = ^j=-dt (2.3) 

Using ([23]) we get E[L{Ci{t))] = VSni. As expected: In = 1. 
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11.2.2 n = {0,7i} 

Now Tq is the exit time of linear Brownian motion from the interval [—1, +1] and 

E[exp(-yTf,)] = l/coshA (2.4) 

with density 

oo ^ 

P(To edt) = n ^(-l)"(n + -)e-("+^)'-'*/2rft (2.5) 

n=0 

and distribution function 

PiTn > t) = - (zl)^g-(n+i)V^t/2 ^2.6) 

^ n=0 ^ ~^ 2 

see e.g. Biane, Pitman, Yor (2001). Using again (12. ip gives 

E[L{C,m = ^vs^t. 

Thus = f • 

11.2.3 fi=[0,27r[ 

The set is the unit circle, and the time Tq is the first hitting time of 1 by a Bessel 
process of dimension 2, therefore 

E(e-4^.)^f| (2.7) 



where 

oo 



°° 2k 



A;=0 

is the Bessel function. Thus 







/o(A)- 

A calculation by Mathematica (thanks to Daoud Bshouty for that) gives = 2,08323.. 



II. 3 A cone of angle 7r/3 

We consider the case Q = {0, 27r/3}. Therefore is the exit time from a cone of angle n/3 
tangent to the unit circle. In this case one can use results from Doumerc and O'Connell 
(2005). There, the exit time from the cone Xi > X2 > in R^, starting from the point 

{xi,X2,X3) is expressed as 

P(r >t)= Pi2 - Pl3 + P23 
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with pij = Jg'*' ^j)/^g y^/'^dy. Taking the orthogonal projection of Brownian motion 
on the hyperplane Xi + X2 + x^, = 0, with the starting point (a/2, 0, — \/2), one gets 

P{T^>t) = ^'-\2 e-y'/'dy - e~y'/'dy 



and a straightforward computation gives 
and 



l/2n _ 3 /2 



Thus £n = 3/2. 



II. 4 An equilateral triangle 

An interesting case is when = {0, 27r/3, 47r/3} and Pn is an equilateral triangle centered 
at 0. We will compute the distribution of and get a curious identity in law. 

In order to do the computation it will be convenient to take the equilateral triangle 
A with vertices 0, 1, — j'^ where j — (—1 + iy/3)/2. Its center is Xq — {1 — i^)/3. The 
case of the triangle Pn can be obtained by an affinity. The orthogonal reflections with 
respect to the three lines bordering the triangle A generate a group W of affine isometrics 
of the Euclidian plane, whose fundamental domain is A. Using the reflection principle, 
we can compute the probability transition of Brownian motion in the triangle, killed at 
the boundary. One obtains 

Pti^^y) = 'iet{w)pt{x,w{y)) 
wew 

where pt{x,y) is the ordinary heat kernel on the Euclidean plane and x and y are points 
of this plane. . If we denote R the lattice generated by 1 and j, then the dual lattice 
R (such that {w,w) G Z) is generated by the dual basis a = 1 + ^ and b = The 
plane is tiled by translates of A by R, which correspond to direct isometrics w &W and 
translates of A. Thus 

Pti^^y) = Ypt{x - y,r) - pt{x - y,r). 

reR 

The function 

Kt^z) = YPt{z,r) 

reR 

is periodic in z with period lattice it! andp°(x, y) = h{t,x — y) —h{t,x — y). One can apply 
Poisson summation formula, i.e. expand it in terms of the characters e^*'^^ '^^ obtaining 

feR 
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Let us now integrate with respect to y in the triangle to obtain the distribution function 
for Tn, the first exit time of Brownian motion starting from xq, 



We use the coordinates y — yi + y2j, 1 > Z/i > 2/2 > on the triangle, with Jacobian 
^ and r = ma + nb on R, then 

{y,f)—myi+ny2, (xq, r) = 2m/3 + n/3, \w\'^ — 4{m'^ + mn + n'^)/3. 

Furthermore, by periodicity 



h{t,x-y)dy= / h{t, x ~ j - y)dy = / h{t,x-y)dy 

J A JA+j 

and the triangle A + j corresponds to the coordinates 1 > 2/2 > yi > 0. Let us compute 

/ hit,Xo-y)dy = ^y2^2^.(2m/3+nm f\yf'\-2Mmy^+ny,)^y^^-SnHm^+mn+n^)t/3 

It is easy to see that terms with riT^O, 777.7^0, m + n^^O give contribution. Also the 
term m — n — Q will be cancelled by the other integral. The sequence m = 0, n 7^ 0, gives 

"1 /•{/! 



The sequence m 7^ 0, n = 0, gives 

_ /■ ^^^^^g-2i.2/img-87r2m2t/3 

= -E 

and the sequence m 7^ 0, m + n = 0, gives 



2i7rm 



E 



-87r2m2t/3 



^2i7rm/3 ^-87r^m^t/3 

2i7rm 
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The sum of these three terms is 

Sm(27rn/3) ^_87r^n^t/3 



-e 

nn 



n=l 

The other integral is 

'1 /•y2 



I h(t, Xo - y)dy = V e2'^(2m/3+n/3) f r ^-2MmyM,)^y^^-8nHm^+mn+n^ 

Using 



)t/3 



Jo JO 



if m or n is 7^ we see that, except for the term m, = n = the other terms are the 
opposite of what we computed, therefore after some rearrangement 



OO / s 



nn 

n=l 



where 

2 

X3{n) = --^sin(27rn/3) 
is the multiplicative Dirichlet character modulo 3. The density is 

CXD 

P(Ta e dt) = 8\/3 5]7rnx3(n)e-«'^'"'*/^ 



n=l 



II. 4.1 A curious identity in law 

Let 5" be the first hitting time of 3a by a three-dimensional Bessel process starting from 
a. The Laplace transform is 

smh 3Aa 

^ ' n=l 

Inverting term by term gives the density 

^g^^_X3(n)e-2"^«V* 



We can also use the alternative expression 



3sinh(Aa) ^> 3\/37rnx3(n) 



sinh(3Aa) ^ n^n^ + 9X^a^ 

^ ' n=l 
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and invert term by term the Laplace transform to get 



18a2 

n=l 

The agreement between these two expressions is an instance of the functional equation of 
theta series. Putting 8/3 = l/(18a^), or a = we see that the exit time Tq and the 
time S have the same distribution. Furthermore, the Mellin transform of this distribution 
is 

E{T^) = 7i-^'8^-'3'+^^^T{s + l)L^3(2s + 1) 

where L^^ is the Dirichlet L-function associated with the character Xs- Coming back to 
Tq, the triangle Pq has sides of size 2a/3 therefore 

^ 3sinh(A/2) 
^ ^ sinh(3A/2) ■ 

Using (12.11) we are reduced to an elementary integral, and the result is 
II. 5 The other cases 

II.5.1 n = {0,n/2} 
We have 

huie) = max{X(s) cos^ + Y{s) sinO, Y{s) cosO - X{s) sinO ; 0<s<t}. 

The trick is to use the fact that {X{s) cos6' + F(s) sin6')s>o and {Y{s) cos9 — X{s) sin^^)s>o 
are two standard one dimensional Brownian motions which are independent. As a conse- 
quence, using (12. 3p we have 



PT{hn{0) <h) = H\h/Vi) 



where 



It leads to 



, , . doc , , 

H{z) = -1 + 2 I e---=. (2.8) 
-oo V 2n 



We write 



/•oo l-OO 

E{L{n))=27i {1- H\h/Vi))dh = 271 Vi {1-H^{z))dz. 
Jo Jo 



H^(z) = 4 / 6 2 ^= 1 - / e " 



27r V Jz V27r 
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By the change of variables x = zu and y = zv and w = z"^ /2 we get easily 

""^ dudv 



(poo fC 
2- J J 



Introduce the function (x, y) i-^ g{x, y) defined for x and ?/ in R \ {0} by 



g{x,y) = ^^^^^. (2.9) 

xy 

A tedious calculation shows that if a < 6 and c < d with (0, 0) ^ [a, b] x [c, d] and if 
furthermore a,b,c,d ^ we have 

(^2 ^ ^2)3/2 = ^(^' ^) + d) - g{a, d) - g{h, c). (2.10) 
From this we get that (-^2^"'2")3/2 = 2 — -\/2 which finally gives that 

in = V2. 

II.5.2 = {0,7r/2,7r,37r/2} 

The method is quite similar to the method used before, but the distribution function H 
of the maximum of a Brownian motion is replaced by the distribution function L of the 
maximum of its absolute value, which can be obtained from (12. 6p . Still with the previous 
notations we have 

/icn(^) = max{|X(s)cose + r(s)sin^|, |r(s) cos^ - X(s) sin^l ; 0<s<t}. 

We use again the fact that (X(s) cos^^ + F(s) sin6')s>o and (Y{s) cos9 — X{s) sin6')s>o are 
two standard one dimensional independent Brownian motions. We get 

Pr{hcM<h) = L^{h/Vi) 

and 

p27r poo roo 

E{L{Cn{t)))= E{hc,it){9))d9 = 2n {1 - L^{h/Vi))dh = 2nVi (1 - L\z))dz. 
Jo Jo Jo 

We now write 

oo oo -(4^+1)2 2 /•(4n+3)2 



1-L''(z) = {l-L{z)){l + L{z)) = - J2 Yl / e-^dx e-"-dy. 

^ k=-oo n=-oo J{ik-l)z i(4n+l)2 



We are led to the computation of the integral for fixed (fc, n) G Z 
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/2 roo I MAk+l)z 2 /•{4n+3)2: 
- / / e~~dx / e'^^dy ) dz. 

^ Jo W(4fc-l)z J{4n+l)z 



By the same change of variables as before x = zu and y = zv and w = z /2 we get 

In order to compute I{k, n) we use the function g defined by (12.91) and the identity (12.101) . 
From ( 12. lip it is apparent that I{k,n) > 0. To see also that 



oo oo 



J2 Hk,n)<oo (2.12) 



— 

k=—oo n=—QO 



we use (12. lip again for observing that in is the integral of the function {u^ + v'^) over 
the union of squares 

D = ur=_oo U^=_^ [4fc - 1, 4fc + 1] X [An + 1, 4n + 3] 

which satisfies D C A = {{u,v) ; u"^ + v'^ > 1}. This implies, using polar coordinates 
u = r cos 6 and v = r sin 6' : 



< / = 27r / r"2rfr = 27r. 



^(^2 + ^2)3/2 (^2 + ^2)3/2 



In order to have some numerical evaluation of £n one can observe that we have the 
following two symmetries: for all (A;, n) G 1? 

I{k, n) = I{-k, n) = I{k, -n-1). 

This implies that Iq can be written as in = + Si with 

So=l^ mn)=4^l ^^^^ l=4L(-l) 

n=—oo n=0 \ / m=0 

(2.13) 

with = ^1 + j2d^~^^ '^i = 4E^=o5]^i^(^)"')- Since {um)m>o is decreasing 

we can estimate ^ ~ 0.3725... We compute by Mathematica (we thank Daoud Bshouty 
for this) and we get 

in = 1,9178... 

II.5.3 Q = {0,n/2,7T} 
Again the same method gives 

/•oo roo 

E{L{Cn))=27r {l-L{h/Vi)H{h/Vi))dh = 2nVi {I - L{z)H{z))dz. 
Jo Jo 

We write for simplicity L = 1 — 2Li and H = 1 — 2Hi so that 

POO 

E{L{Cn)) = 2TTVi / (2Li(z) + 2Hi{z) - ALi{z)Hi{z)) dz. 
Jo 
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A previous calculation has shown that 2TX\ft 2Li{z)dz = ^y/Sirt and Letac (1978) -or 
a direct calculation- shows 2-K\/t 2Hi{z)dz = VSni. The delicate remaining point is 
the computation of 4:Li{z)Hi{z)dz. By the same methods as before, the same changes 
of variable x = zu and y = zv and w = z^/2, the use of the function g of (12.9p and the 
formula (I2.10p we get 

/•oo °° POO p 

27iy/t / ALi{z)Hi{z)dz = VS^ / / 



n=— oo 
oo 



°° dudv 

1 .4n+l W+W^ 



= Vs^t J2 (9{h 4n + 1) - g{l, An + 3)). 

n=— oo 

It is easily seen that with the notation (I2.13P 

J] {g{l, 4n + l)- g{l, An + 3)) = ^ ^ 0.7450... (2.14) 



- + 0.2550. 



Thus 

in 

z 

III Comments 

1. The explicit calculations of the paper for the mean perimeter are feasible since we are 
able to have the exact distribution of the support functions 9 i— > hc{0) for the random 
convex sets considered here. The general case seems a difficult problem. 

2. Calculations for the mean areas rely on the formula A{C) = /q (/ic(^)^ - h'c{ef)de, 
where h'(j{6) denotes the derivative on the left, which can be shown to exist (see Letac 
(1983)) and ultimately of the distribution of h'(j{6) when C is random. This problem was 
solved in the unpublished thesis of El Bachir (1983) for Ci(t), where E(A(Ci(t)) = 7rt/2 
is proved. It is unsolved for the other convex sets considered here, except for Cooii) the 
smallest circle, centered at 0, surrounding Ci(t), indeed 

ds 

¥.{A{Coo{t)) = 7Tt / =- = vrt X 3.06883. 

Jo /o(v2s) 

a little more than six times the quantity K{A{Ci(t)). Thanks to Daoud Bshouty for the 
numerical evaluation of this integral. 
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